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We analyze the spectral properties of a very general two-channel fermion-boson transport model 
in the insulating and metallic regimes, and the signatures of the metal-insulator quantum phase 
transition in between. To this end we determine the single particle spectral function related to 
angle-resolved photoemission spectroscopy, the momentum distribution function, the Drude weight 
and the optical response by means of a dynamical (pseudo-site) density-matrix renormalization group 
technique for the one-dimensional half-filled band case. We show how the interplay of correlations 
and fluctuations in the background medium controls the charge dynamics of the system, which is a 
fundamental problem in a great variety of advanced materials. 
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I. INTRODUCTION 

Charge transport normally takes place in some back- 
ground medium. To understand how the environment 
affects the moving carrier and vice versa is a difficult 
question and in this generality at present perhaps one 
of the most heavily debated issues in condensed mat- 
ter physics. Here the term "background" describes a 
variety of situations. We can think of the motion of a 
hole through an ordered insulator-i. Examples are the 
high-Tc cuprates and the colossal magnetoresistive man- 
ganates, with a background of spins and orbitals, respec- 
tively, forming a pattern of alternating order. Then, as 
the hole moves, it disrupts the order of the background, 
which on its part hinders the particle transfer. Never- 
theless coherent particle transport may occur, but on a 
strongly renormalized energy scale. The new quasiparti- 
cles formed in the cuprates and manganates are spin or 
orbital polarons^i^. Another situation concerns a charge 
carrier coupled to a deformable background. Here, if the 
interaction with phonons is strong, the particle has to 
carry a phonon cloud through the medium. The out- 
come might be a "self-trapped" small lattice polaron^. 
In this case hopping transport, accompanied by phonon 
emission and absorption processes, evolves as the domi- 
nant transport channel. 

So far we have considered a single particle only. It 
is quite obvious that the problem becomes even more 
involved if the particle density increases. Then the inter- 
relation between charge carriers and background medium 
may drive quantum phase transitions. The appear- 
ance of ferromagnetism in the three-dimensional (3D) 
manganates, superconductivity in the quasi-2D cuprates, 
or charge-density-wave (CDW) states in ID halogen- 
bridged transition metal complexes are prominent exam- 
ples^. In the theoretical description of these strongly 
correlated systems an additional difficulty arises: The 
particles which are responsible for charge transport and 
the order phenomena of the background are the same. 
As a consequence, on a microscopic level, rather involved 
many-particle models result, which incorporate the cou- 



pling between charge, spin, orbital and lattice degrees of 
freedoroiii. Naturally this prevents an exact solution of 
the problem even in reduced dimensions. 

II. MODEL AND METHOD 

A way out might be the construction of simplified 
transport models, which capture the basic mechanisms 
of quantum transport in a background medium in an ef- 
fective way. Along this line a novel quantum transport 
model has been proposed recently^'^,^ 

^ =-^'' E fUM+b,)-^T.^bi+b^)^oJ2^h , (1) 

which mimics the correlations inherent to a spinful 
fermionic many-particle systems by a boson affected hop- 
ping of spinless particles oc tb (see Fig. [1]). In the 
model (P), a fermion f^^^ creates (or absorbs) a local 

boson bf"^ every time it hops, which corresponds to a 
local excitation in the background with a certain en- 
ergy lvq. Because of quantum fluctuations the distor- 
tions are able to relax oc A. A unitary transformation 
bi hi + X/loq replaces this term by second trans- 
port channel Hf = ~tfJ2{i j) fjfi^ describing unaffected 
fermionic transfer, however with a renormalized ampli- 
tude tf ^ 2\t\j/uiQ. It has been shown^ that coherent 
propagation of a fermion is possible even in the limit 
A = = by means of a six-step vacuum- restoring hop- 
ping process, 

^f' = -^I+2-^I+l-^l^*+2^i+l-Ri ; (2) 

where r\ - flh+A and l\ - fU,_A- Note that 
i?!^"* acts as direct next neareast neighbor (NNN) transfer 
"fi+ifi" ' in complete analogy to the "Trugman path" of 
a hole in a 2D Neel-ordered spin background^". 

The model ([T]) has been solved in the single particle 
sector (A^e — 1) by exact diagonalization^, using a ba- 
sis construction for the fermion-boson (many-particle) 
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Hilbert space that is variational for an infinite lattice 
(A^ = (X))ii. The transport behavior was found to be sur- 
prisingly complex, reflecting the properties of both spin 
and lattice polarons in t-J- and Holstein-type models. 

For the ID half-filled band sector (TVg = iV/2), evi- 
dence for a metal insulator transition comes from small 
cluster diagonalization&lS,. Quite recently the ground- 
state phase diagram of the model ([T]) has been mapped 
out in the whole A-wq plane^-, using a density-matrix 
renormalization group (DMRG) technique^"*. A quan- 
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FIG. 1: (Color online) Schematic diagram of quantum trans- 
port in a background medium. The background could repre- 
sent a magnetically, orbitally or cliarge ordered lattice, but 
also a heat bath or certain chemical side groups. Then the 
proposed transport model H} describes a very general situa- 
tion: As a charge carrier (•) moves along a ID transport path 
it creates an excitation with energy luq (*) in the background 
medium at the site it leaves or annihilates an existing exci- 
tation at the site it enters. It is a plausible assumption that 
the (de)excitation of the background can be parameterized 
as a bosonic degree of freedom . In the case of spin devi- 
ations, orbital fluctuations or lattice vibrations, the bosons 
might be viewed as a Schwinger-bosons, orbitons or phonons. 
Of course, any distortion of the background can heal out 
by quantum fluctuations. Accordingly the A-term allows for 
spontaneous boson creation and annihilation processes. The 
upper panel displays the single-particle case. Depending on 
the model parameters quasi-free, diffusive or boson-assisted 
transport takes placed. The latter case corresponds e.g. to 
the motion of a hole through an ordered antiferromagnetic in- 
sulator. The lower panel shows the half-filled band case. Here, 
for spinless fermions in ID, a repulsive Tomonaga-Luttinger 
liquid evolves, provided the excitations of the background are 
energetically inexpensive (tJo < wo,c) or will readily relax 
(A > Ac(ii'o)). This defines the fluctuation dominated regime. 
By contrast, strong background correlations, which develop 
for large ujo and small A <C it tend to immobilize the charge 
carriers and even may drive a metal insulator transition by 
establishing CDW long-range orde r^^'^^ . 



tum phase transition between a Tomonaga-Luttinger liq- 
uid (TLL) and CDW was proven to exist. A complemen- 
tary study of the dynamical properties of the system is 
therefore desirable. 

In the present work, we employ the dynamical DMRG 
(DDMRG) methodic in order to investigate the effects 
of background fluctuations and correlations on the dy- 
namics of charge carriers in the framework of the ID 
half-flUed fermion-boson model IJ). Thereby the fo- 
cus is on the wave-vector resolves single-particle spec- 
tral function probed by angle-resolved photoemission 
spectroscopy (ARPES) and on the optical conductivity 
probed e.g. by reflectivity measurements. 

In general the dynamic response of a quantum system 
described by a time-independent Hamiltonian H is given 
by the imaginary part of correlation functions of type 



^oH= lim-(^o|Ot — 

1)^0 TT (£/o 
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rOIV'o), (3) 



where the operator O identifies the physical quantity of 
interest. and Eq give the ground-state wave function 
and energy of H. The small 77 > shifts the poles of the 
related Green function Go(w-|-ir?) into the complex plane. 

Single particle excitations associated with the injection 
or emission of an electron with wave vector fc, A+(fc,a;) 
or A^{k,uj), can be written in the spectral form 



(4) 



where = fl, ff7 = f^- iV'o) is the groundstate of a 
iV-site system in the A"e-particle sector while 1-0^) denote 
the n-th excited states in the ± 1-particle sectors with 



excitation energies = — Eq. 

Optical excitations, on the other hand, connect states 
in the same particle sector with a site-parity change. For 
a system with open boundary conditions (OBC) the reg- 
ular part of the optical absorption 

(Treg(tj) = ^^u;n\{iJn\P\lpo)\'^ S[UJ - UJn] (5) 



is related to the dynamical polarizability, (Treg(w) = 
uja{u}), where P = — 3 ifjfj ~ ^) ^he dipole op- 

erator (in units of e), and tOn — {En — Eq). Then the 
current operator is obtained from J = i[H,P]. Applying 
periodic boundary conditions (PBC), the optical conduc- 
tivity can be calculated from 



^ reg 



KV-nlJl^o)!' 



(6) 



Note that for our fermion-boson model ([T]), the current 
operator has two contributions, J = J/ + Jf,, where 
Jf - ^tf flj, - and J, = it, flj, b] - 



/J/;+i^- + fLifjb] - fU,-i^- The f-sum rule 



Sreg(oo) 



ttD 



-TTEktn/2 



(7) 
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connects the frequency-integrated optical response 



5, 



regy^i — ^reg{^') d^' to the kinetic energy i?fei„ = 

;^(0|i? — woX^i where the Drude part oc D serves 

as a measure for coherent transport. For OBC, only a D 
precursor exists in the metallic region. 

In the actual DDMRG calculation of spectral functions 
the required CPU time increases rapidly with the number 
of the density-matrix eigenstates m. Since the DDMRG 
approach is based on a variational principlei^, we first 
of all have to prepare a good "trial function" for the 
ground state with as many density-matrix eigenstates as 
possible. As a rule we keep m ~ 500 states to obtain 
the true ground state in the first five DDMRG sweeps 
and afterwards take m ^ 200 states for the calculation 
of the various spectra from ([3]) with a broadening 77 = 
0.1. In order to save CPU time in the DDMRG runs we 
take into account just nf, = 3 pseudosites. In this case 
the rifj-th local boson pseudo-site density is smaller than 
10~^. Using nf, = 4 this value can be reduced to 10"^ 
which leads, however, not to visible change of the spectra 
because the discarded weight in the DDMRG calculations 
is ~ 10~^ (i.e. three orders of magnitude larger than for 
the DMRG ground-state calculations). 



III. RESULTS 
A. Photoemission spectrum 

Let us first discuss the single-particle spectra of the 
transport model (H)) in the regime where the background 
is stiff, i.e. the distortions induced by the particle hop- 
ping process are energetically costly (wq = 2). 

For very large A the free transport channel nevertheless 
dominates and an almost particle-hole symmetric spec- 
trum (yl+(fc,tt' — Ep) ~ A~{k — TT, Ep — ijj)) results (see 
Fig.[21upper panels). As A decreases, the background dis- 
tortions hardly relax. Consequently the bosonic degrees 
of freedom will strongly affect the transport. The middle 
panels of Fig.[2]show how, at A = 0.1, strong correlations 
develop in the occupied states probed by photoemission 
(PE) for w < Ep- The introduced hole can only move 
coherently by the 6-step process where in step 1-3 
three bosons were excited, which are consumed in steps 
4-6 afterwards. In this way the collective particle-boson 
dynamics leads to a flattening of the "coherent" band 
for k < kp. By contrast an additional electron, which 
probes the unoccupied states in an inverse (I)PE experi- 
ment (w > Ep), can more easily move by a two-step pro- 
cess, even if pronounced CDW correlations exist in the 
background mediumi^. The incoherent parts of A{k,uj) 
far away from the Fermi energy Ep are caused by ex- 
citations with additional bosons involved (bear in mind 
that the ground state with Ne electrons is a multi-phonon 
state and the wave vector of the Ne_± \ target state corre- 
sponds to the total momentum of electrons and bosons) . 

While for A — 0.1, A{kp,ijj) has finite spectral weight 
at Ep, i.e. the system is still metallic (albeit the TLL 
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FIG. 2: Intensity (left panels) and line-shape (right panels) 
of the single-particle spectral function A{k,uj) in the half- 
filled band sector of the fermion-boson transport model ((TJ 
on a A' = 32-site chain. The upper two rows (lower row) 
give DDMRG data for A = 2 and 0.1 in the metallic regime 
(A = 0.01, insulating regime), where wo = 2. All energies 
are measured in units of it,. Since we apply OBC, we use 
quasi-momenta k — ■nl/{N -I- 1) with integers 1 < / < A. 



charge exponent Kp is noticeably reduced from onei^) , an 
excitation gap opens in the PE spectrum as A falls below 
a certain critical value, provided that loq > a;o,c(A = 0)^'^. 
We find Xc{coo = 2) ~ 0.05. The lower panels of Fig. [H 
show A(k, Lu) for A = 0.01, in the insulating regime, where 
a CDW with true long-range order exists. The TLL- 
CDW quantum phase transition is driven by the corre- 
lations that might evolve in the background medium at 
commensurate fillings. Let us emphasize the dynami- 
cal aspect of this process: The (collective) bosonic ex- 
citations are intimately connected to the motion of the 
particles, and themselves have to persist long enough in 
order to affect the many-particle state. The ARPES 
spectrum for the insulating state clearly shows the dou- 
bling of the Brillouin zone. The remaining asymmetry 
with regard to the spectral weight of the absorption sig- 




io — Ep uj — Ep 

FIG. 3: Density (left) and line shape (right) plot of the A{k, u) 
spectra. Again TV = 32 (OBC), A = 0.01, but now luq = 1- 

nals as fc ^ (tt — /c) vanishes for A — > 0. Most notably 
the widths of the highest PE and lowest IPE coherent 
bands differ by a factor of about (ib/wo)'*, since the 
CDW order is restored if the injected hole [electron] is 
transfered to a NNN site by a process of order 0(<^/wq) 
[0{tl/ujQ)]. Hence the CDW state exhibits a correlation- 
induced asymmetric band structure^^. 

The strong interrelation of charge dynamics and back- 
ground fluctuations becomes obvious if we decrease ujq 
below wo.c keeping A = 0.01 fixed. Of course, in passing 
the accompanied insulator-metal transition the PE spec- 
trum changes completely, but the "nature" of the TLL at 
Wo = 1 is different compared to that of the metallic state 
realized at larger ujq and A as well (cf. Fig. [3] and upper 
panels of Fig. [2]). The single-particle spectrum for luq = 1 
shows sharp absorption signals in the vicinity of kp only. 
In a wide fc-space region emanating from fc = (fc = tt) 
the PE (IPE) spectrum is smeared out (over-damped), 
i.e. here the dynamics of the system is dominated by 
bosonic fluctuations. 



B. Momentum distribution function 

The different transport behavior becomes also appar- 
ent in the momentum distribution function 

By means of DMRG, the ground-state correlation func- 
tion (CjCi) can be easily calculated for PBC. Figure [4] 
displays n{k) for two characteristic boson energies, above 
and below loq^c- 

In the former case, the TLL-CDW transition causes 
significant changes in the functional form of n(k). For 
A > Ac, one expects an essential power law singularity 
at kp^ corresponding to a vanishing quasiparticle weight. 
For finite TLL systems the difference A = n{kp — S) — 
n{kp + 5) is finite (with 5 = 7r/66 in our case)^^. A 
rapidly decreases approaching the CDW transition point 
with decreasing A (see data for ujq = 2 (red squares)). In 
the CDW phase the singularity at kp vanishes. Note that 
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FIG. 4: (Color online) Momentum distribution function n(fe) 
for the half-filled two-channel transport model ([T| with 66 
sites and PBC, as A decreases from one (a) to zero (d) at 
cjo = 1 (circles) and 2 (squares). Triangles show n(fc) for a 
half-filled t-t'-U model (38 sites, PBC) with U = lOt', t = t' 
[panel (c)] and t = Q [panel (d)] (see text). In order to obtain 
more accurate ground-state data we use nt = 4 pseudo-sites. 

the periodicity of n{k) doubles at A = 0, in accordance 
with a NNN-only hopping channel. To substanti- 
ate this reasoning we have included in Fig. 0] 7i(fc)~data 
calculated for the ID Hubbard model with additional 
NNN-transfer t' . We see that n{k) of the fcrmion-boson 
model (U)) is in qualitative agreement with our data and 
previous results for the t-t' Hubbard model^^, in partic- 
ular for the case t = Q. The upturn in n{k) for k > kp 
persists even in the metallic regime as long as NNN hop- 
ping processes triggered by the (CDW) correlations in 
the background are of importance. 

For Wo = 1 the system stays metallic for all A. Be- 
sides the usual renormalization of n{k) with increasing 
correlations (i.e. decreasing A) we find a slight upturn in 
n{k) for k < kf. This might be attributed to the fact 
that in our model ([1]) a particle injected with k — ±7r is 
almost unaffected by bosonic fluctuations (which holds 
also for the single-particle case^). So to speak the sys- 
tem behaves as a nearly perfect metal at this point. It 
is worth to mention that an increase in n{k) for both 
k "^kp and k > kp has been found also for the momen- 
tum distribution function of the Hubbard model (with 
and without magnetization) using the Gutzwiller varia- 
tional wave functioni^. 



C. Optical response 

Finally we consider the evolution of the optical conduc- 
tivity going from the correlated TLL to the CDW phase 
at cjQ = 2. The corresponding optical absorption spectra 
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FIG. 5: (Color online) Optical conductivity areg{uj) and in- 
tegrated spectral weight Sreg{i-o) in the half-filled transport 
model lU. The main panel (inset) displays DDMRG data for 
ti^o = 2 at A < Ac (A 2> Ac), i.e. in the insulating (metallic) 
region, applying OBC (PBC). 

are depicted in Fig. [5] In the metallic state most of the 
spectral weight resides in the coherent Drude part. At 
A = 1 (see inset), we find ttD/N ~ 1.6, which has to 
be compared with Sreg{oo) — 0.2 (of course D decreases 
as A gets smaller). In this case the wave- vector resolved 
single-particle spectra roughly extends from a; = — 6 to 
w = 6. The regular part of the conductivity is mainly 
due to excitations to the phononic side bands appearing 
in the sectors with momenta far away from kp. In the in- 
sulating region (see main panel), the first peak at about 
uj ~ 0.5 can be assigned to an optical excitation across 
the gap in the (coherent) two-band structure. These ex- 
citations are only accessible for A > 0. Additional ex- 
citations with higher energy occur around multiples of 
the boson frequency, where w ~ = 2 sets a absorp- 
tion threshold for the A = case. As expected for an 



insulating system with OBC, the whole spectral weight 
is contained in S'reg(oo) ~ —irEkm/'^- We emphasize 
that the CDW state in our model contains less than one 
boson per site on average, unlike e.g. the Peierls insu- 
lating state in the Holstein model. That is the CDW 
phase typifies rather as a correlated insulator — such as 
the Mott-Hubbard insulator — and no divergence occurs 
at the optical absorption thresholdi^. 



IV. SUMMARY 

In conclusion, we have determined the spectral prop- 
erties of a highly non-trivial two-channel fermion-boson 
transport model for the ID half-filled band case, using an 
unbiased DDMRG technique. The background medium, 
parameterized by bosonic degrees of freedom, strongly 
influences the charge carrier dynamics, as it happens in 
many novel materials. If the background fluctuations 
dominate we find diffusive transport. In opposite case of 
strong background correlations coherent quantum trans- 
port may evolve on a reduced energy scale. These corre- 
lations can also trigger a metal insulator transition. The 
insulating CDW state has a asymmetric band structure, 
leading to characteristic signatures in the ARPES and 
optical response. Whether an extended model with spin- 
ful fermions gives rise to an attractive metallic phase like 
in the Holstein-Hubbard model^° would be an interesting 
question for further research. 
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